Raman scattering on phonon— plasmon coupled modes in magnetic fields 
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Raman scattering on phonon-plasmon coupled modes in high magnetic fields is considered theo- 
retically. The calculations of the dielectric function were performed in the long-wave approximation 
for the semiclassical and ultra-quantum magnetic fields taking into account the electron damping 
and intrinsic lifetime of optical phonons. The Raman scattering has resonances at the frequencies 
of coupled modes as well as at multiples of the cyclotron frequency. The dependence of the Raman 
cross section on the carrier concentration is analyzed. 



I. INTRODUCTION 

Phonon-plasmon coupled modes arise as a result of 
the interaction of optical phonons with charge carriers 
in semiconductors. These modes are usually observed in 
the absence of a magnetic field, for instance, in Raman 
scattering Q . Interest in these modes is related not only 
to applications in optics but also to principal questions of 
condensed-state physics, in particular, to the possibility 
of determining the magnitude of electron-phonon inter- 
action with their help. The first experimental results for 
the Raman scattering on phonon-plasmon coupled modes 
in magnetic fields were recently obtained ^ , but only the 
preliminary theoretical paper 0] has been published up 
to now. 

The longitudinal optical phonons interact strongly 
with free carriers because of an electric field accompanied 
in the lattice vibrations. Let us estimate the parame- 
ters of the corresponding electron-phonon system. The 
frequency of an emitted or absorbed optical phonon u is 
usually equal to several hundreds of degrees, for example, 
this is 35 mV in GaAs. If the light scattering is excited by 
the laser frequency uji — 1.5 eV, the momentum transfer 
k from the incident radiation has the order of Ui/c ~ 10^ 
cm^^. Then, the conditions uj/c <C k <C oj/vp are ful- 
filled, where the Fermi velocity is vp ~ 0.5 x 10* cm/c 
in the case when the statistics of carriers is degenerate. 
The left-hand side of the above inequality means that 
the electric field associated with the phonon vibration 
of frequency uj is really static and therefore the electric 
field is longitudinal with respect to the wavevector k. 
The right-hand side of the inequality allows the electron 
contribution to the dielectric function £(k, u) to be calcu- 
lated using a series expansion in terms of the dispersion 
parameter kv/uj. We shall see that at least the lowest- 
order correction should be held, because it has a resonant 
character. As for the lattice contribution to e(k, w), it 
can be taken at /c = 0, since the phonon dispersion is 
negligible for such the momentum. 

The effect of magnetic field on optical phonons is most 
essential if the cyclotron frequency of carriers ojc has the 
order of the optical phonon frequency. This means that 
the magnetic field must reach 20 T for effective carrier 
mass of 0.063 mo in the same GaAs. Changing the car- 
rier concentration, we can observe both the semiclassi- 
cal regime, when the cyclotron frequency is much less 



than the Fermi energy [ojc <ti sp), and the ultra-quantum 
regime, when ujc > £f- 

Two circumstances should be noted. Firstly, both the 
frequency of the coupled modes and the width of the cor- 
responding resonance are of interest in the Raman light 
scattering. Therefore, the contribution of carriers to the 
dielectric susceptibility must be calculated with regard 
to their damping. The spatial dispersion of susceptibility 
in a magnetic field under these conditions has not been 
calculated so far. Secondly, in addition to the Coulomb 
interaction of carriers with the longitudinal phonon vi- 
brations, which is taken into account by the dielectric 
function, there exists a deformation interaction with both 
the LO and TO modes, which bears the Frohlich name in 
theoretical works. It arises because of the nonadiabatic- 
ity of electron-phonon systems and leads to a certain (as 
small as the nonadiabaticity parameter) renormalization 
of phonon frequencies. In the absence of a magnetic field, 
this renormalization has been recently considered in - 
6], and it will not be taken into account here. 

The structure of the paper is the following. A short 
outline of the theory of Raman scattering on the phonon- 
plasmon coupled mode is give in Sec. II. Then the di- 
electric function of the system in semiclassical and ultra- 
quantum regimes is evaluated in Sees. Ill and IV. The 
theoretical Raman spectra are presented for the various 
magnetic fields and carrier concentrations in Sec. V. Fi- 
nally, the conclusions are summarized in Sec. VI. 



II. INELASTIC LIGHT SCATTERING ON 
PHONON-PLASMON COUPLED MODES 

Let us consider the Raman scattering on optical vibra- 
tions in a polar lattice with free carriers. We use the no- 
tations bj for the phonon displacements of the branch j. 
The subscript j denotes the various phonon modes: longi- 
tudinal or transverse ones. More precisely, the subscript 
j indicates the different phonon representations which 
can be degenerate. The transformation properties of the 
coupling constants gj are determined by this representa- 
tion. The LO vibrations are accompanied by the internal 
electric field E(r,t). We consider the phonon displace- 
ments and electric field as classical variables because the 
magnetic field affects quantum-mechanically only the free 
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The effective Hamiltonian describing the inelastic hght 
scattering on phonon-plasmon modes can be written as 



H 



cPrN{v, t)U{r, t), 



where 



M{r,t)=gjb,{r,t)+gEE{r,t) 



(1) 



(2) 



is a hnear form of variables bj and E. We do not in- 
clude the pure electronic Raman scattering by free carri- 
ers, which has been considered in many papers (see, e.g. 
0,0). 

The notation U{r, t) is introduced for a product of the 
vector-potentials of incident and scattered photons: 

A^'\r,t)A^'\r,t) = C/(r, i) = exp[z(k • r - tjt)]t/(k, w), 

where the momentum and frequency transfers k = k*-*-* — 



k(^), 



UJ 



UJ 



The polarization vectors e'a and 



ep of the incident and scattered photons are included in 
the coupling constants as well as the polarization vector 
of E(r,t) is included in the coupling constant gs- The 
estimation of the coupling constants gives gj ^ /a'^, gs ^ 
1/ea where a is the lattice parameter. 

If we evaluate the generalized susceptibility x(k, uj) in 
the linear response to the force U (k, ui) 



A/-(k, 



^x(k,c^)C/(k,^), 



(3) 



we obtain the Raman cross-section 



da 



7rc(l - exp(-w/T) \chmu!^'-\ 

x|C/(k,w)pImx(k,w), (4) 

(s) 

where kz is the normal to the sample surface component 
of the scattered wave vector in vacuum. 

One note should be made here. Of course, any sample 
has the surface. The surface effects in the Raman scatter- 
ing was considered in the work and they are omitted 
in the derivation of Eq. Furthermore, the incident 
and scattered fields do not penetrate into the bulk due 
to the skin effect. For the optical range of the incident 
light, we have the normal skin-effect conditions. Then 
we integrate in Eq. Q the distribution |/7(k, over 
the normal component kz- As shown in the paper the 
integration of |[/(k, gives a factor I/C2, where (2 is 
expressed in terms of the wave- vector components inside 
semiconductor: C2 = Im(fci*-' -I- fci**^). The Raman cross 
section (0J obtained is dimensionless. It represents a ra- 
tio of the inelastic scattered light energy to the incident 
energy. 

The equation of motion for the LO mode 6lo has the 
form 



M 



gTot/(k,cj) 
M'N 



(5) 



where N is the number of unit cells in 1 cm'^, M' is 
the reduced mass of the unit cell, and Z is the effective 
ionic charge. Notice that the optical p honons always have 
the so-called natural width F ~ ujQy^m/M . The natu- 
ral width results from decay processes into two or more 
acoustic and optical phonons. In the final expressions, 
we will substitute lu^q — uj"^ ^to ~ ~ '^'^ ■ '^^^ 
equation |SJ| is applied as well to the transverse phonons, 
but the electric field has to be neglected in the case of 
the TO phonons. 

The electric field E(r, i) can be obtained from the Pois- 
son equation divD = 0. There are several contribu- 
tions in the induction _D: (1) the polarization aE{r,t) 
of the filled electron bands, (2) the lattice polarization 
NZbi^o{^,t), (3) the contribution of free carrier density 
p = -divPe, and (4) the term P = -dU/dE = -gEU 
explicitly results from the Hamiltonian, Eqs. JQl, l(2Jl. 
Collecting all these terms into the Poisson equation, we 
find 



eoo£^(k,w)+47riVZ6Lo(k,c^) 

AlTl P 

-^p(k,c.)-47rg£;7(k,c.)=0, (6) 



where the high-frequency dielectric constant = 1-1- 
47ra. 

We shall calculate the carrier density p(k, lo) in the 
following sections and find the electron contribution in 
the dielectric function £e(k, w) = £00 + 47riep(k, 
Then Eq. © takes the fo rm 

£e-B(k, uj) 47riVZ6Lo(k, oj) - 47r5£;[/(k, uj) = 0. (7) 

Solving Eqs. Q, and using Eq. (0), we find A/'(k,a;) 
and obtain the susceptibility defined by Eq. (Q: 



X(k, uj) = Angl 



£e(k, Lj)C'^uj^Q/eooUj'^, - A - 2Cljto 

£e(k,w)A + £ooC^pj 



(8) 

where we set A = "-^to ~ ~ iloT, the plasma frequency 
of ions = inNZ'^ /soaM' , and the Faust-Henry coef- 
ficient C = gToZ/gsM'tj^Q. 

We see from Eq. (|HJ| , that the poles of the general sus- 
ceptibility coincide with the zeros of the total dielectric 
function 



£(k, t^;) = £e(k, cj) 



(9) 



Therefore, the peaks of the Raman cross section give the 
frequencies of coupled modes determined by the condi- 
tion £(k, Lj) — 0. 

In order to calculate the generalized susceptibility, Eq. 
©, we must find the dielectric function £e(k, w), i.e. the 
free carrier density p(k,uj). 
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III. THE DIELECTRIC FUNCTION IN 
SEMICLASSICAL REGIME 

Consider the geometry when the magnetic-field efi^ect 
is most important: the magnetic field H is directed along 
the axis z, whereas the wavevector k and the correspond- 
ing electric field E are directed in the perpendicular di- 
rection X. In semiclassical conditions Uc ^ £f, we can 
use the Boltzmann equation 



i{uj - kvx)f + w, 



dip 



if- < f >)/r, (10) 



written in the t— approximation within the momentum 
variables: the component pz along the magnetic field, the 
angle (p of rotation around the magnetic field, and the 
energy e. The element of phase volume in these variables 
is mdsdpzdip. For a quadratic electron spectrum, we can 
assume, for example, that Vx = v±{pz) simp. The factor 
eEi5(e — e^), containing the charge, the electric field, and 
the Dirac function, is separated out of the distribution 
function /. Angle brackets designate the average over 
the Fermi surface 



{...)mdpzd(p) / I mdpzdip 



< ... >= 



The term with averaging in the Boltzmann equation is 
necessary for the fulfillment of the conservation law of 
the electric charge. 

We assume kv/uj <C 1 and expand the solution to Eq. 
(|10|l in powers of this parameter to the second order / = 
/o + /i + /2. Because (vx) vanishes in integrating over 
(p, it is seen from Eq. H1U|) that (/o) — 0. Therefore, the 
system of equations for /o, /i, and /2 takes the form 

- ILO /o +'^c-r- = Vx, 

dip 

dp T 

where j = 1, 2 and lo* = lo + i/t. The zeroth- and first- 
order approximation, which have to be periodic in the 
angle p, are readily found 

fa= Vx exp [-iu!*{p} ~ p>')/i^c] = 



(11) 



(< /i > /t) - ikvxfo) exp [-iu!*{p - p')/u!c], 

(12) 

First, (/i) and, then, /i can be found by averaging both 
the sides of the last equation. We obtain 



/: 



-ikv , 



Slot 



+ 



A{uj1 - uj* 
{u:c+l^*)v\ 2iv , {uJc-uj*)v\ 



(2wc — uj*)v'l 



-2iip 



(13) 



For the second-order approximation, an expression 
similar to /i, Eq. H12II is obtained. The difference is 
connected with the fact that /2 is an odd function of ve- 
locity, and, therefore, the mean value (/2) vanishes. It 
can be found that 

k 2^ 



/2- 



k^vj_Vp 



8(C.2 



+ 



\3lut 

{uJc+UJ*)v1 
(^c - 



— LO* ^ 
UJ*)v 



UJc — iiJ 



J.ip 



3wc 



Wc — W 



1(14) 



{2ujc + uJ*)v 

It is simpler to calculate the current j(k, cj) instead of 
the electric charge, using the conservation law p(k, lu) = 
kj(\i,Lu)/uj. The contribution of carrier into the conduc- 
tivity is 

2dipdp 



Vxfm- 



(15) 



(27r?i)3 ■ 

Because the first-order approximation /i, Eq. (|13|l is 
even with respect to velocity, it makes no contribution 
to the conductivity. Using the zeroth-order, Eq. I|ll|) 
and the second-order, Eq. H14() approximations, we find 
the carrier conductivity, Eq. (|15|l and the susceptibility 
Airia/uj. With regard to the contribution of filled bands, 
the electron susceptibility can be written as 



,{uj,k) 



Uj{ljJ*'^ — u"^) 



5i 

9lut 



(16) 



The expression obtained remains valid both in the ab- 
sence of a magnetic field [lOc = 0) and in the coUisionless 
limit (r = oo). If both these conditions are fulfilled, Eq. 
(fT?)|l converts to the known expression with the true co- 
efficient 3/5. This coefficient, as well as the others in Eq. 
(|16|l . was calculated here for the quadratic electron spec- 
trum; however, the dependence itself on the frequency 
w, magnetic field, and damping is retained in the gen- 
eral case. It is easy to see that the highest resonances 
nWc contribute to the dielectric function the terms of the 
order (kvp/i^c)"' with the even n. 

The dielectric function of a nondegenerate electron 
plasma was obtained in Ref. |10| 



e{uj,k) = 1 



1- A 



A 



(17) 



where A = {kvth/^jJcf' : and vth is the thermal velocity of 
electrons. At r = cx) and A = [kvp / i^cY it coincides 
with the electronic term in Eq. H16|l. 



IV. ULTRA-QUANTUM MAGNETIC FIELDS 

Now we consider the very high magnetic field when 
only one lowest electron level is occupied. In order to 



4 



calculate the dielectric function, we use the expansion of 
the electron density 



(18) 



in terms of the eigenfunctions 

where s is the total set of the quantum numbers s — 
{n,py,pz} and Xn can be expressed in terms of the Her- 
mite polynomials. The magnetic field is chosen in the 
z— direction. 

Using the equation of motion for the operator afa^i in 
the weak field ^(r), we find the linear response 



e4>s'Afs - fs') 



(19) 



where 4)s' ,s is the matrix element of the potential (/'(r), fs 
is the Fermi distribution function depending on Es and 
the average is taken over the ground state. Inserting 
Eq. p9|l into Eq. p8|l . averaging and taking the Fourier 
transform with respect r, we obtain the induced charge 
and then the electronic dielectric function 



ee(k, w) 



47i-e2 ^ I dxdpydpz ^^k^^ 



fc2 



n,n' 



(27r?i)2 



xXni-cpy/eH)xnix ~ cpy/eH) (20) 

ifs-fs') 



>^Xn'{-cp'/eH)xn'{x - cp' JeH) 



huj — £s' + Es 



where p'^ ^ py + ky, Ss = [n + 1/2 ± g/4)hwc +pl/2m*, 
e^i = (n' + 1/2 ± g/A)huJc +pl/'2m*. In the following, we 
assume that the g-factor equals 2 as for the free electron. 

Let us introduce the new variables: x — cpy/eH 
X, —cpy/eH x' . Then the integral 120|l takes the form: 



dxdx'dp 



xXn{x)xn{x')xn'{x - cky / eH)xn' {x' " cky / cH) 

„ ifs-fs') 



(21) 



Because H is perpendicular to k, we can choose ky = 
kz = 0. Then we get under the integral (|21|l the matrix 
element = {e^''^^)n,n' squared. 

The wave function Xnix) changes over the magnetic 
length qh — y/Ti/m*uJc- We have k ^ l/a_f/ for the 
high magnetic fields. Therefore we can evaluate the in- 
tegral In,n' expanding the exponent in powers of kx and 
using the known matrix elements x„.„_i = Xn-i.n = 
\J nTi/2m*ujc- Due to the Fermi factors fs,fs', the inte- 
grand in Eq. I|21|l does not vanish only if n n' and one 
of the numbers n, n' have to be zero in the ultra-quantum 
limit. Then we are interested in the off-diagonal matrix 



elements /„.„'. Under these conditions, we obtain to the 
forth order in k: 



n',n+l 



\2to*w, 
1 



-{n + l){2n + 3)6n,,n+i 



(n + l)(n-F2)<5„.,„+2 



+ {n^n'). (22) 



Substituting Eq. into Eq. l(?T)) , taking terms with 
n = or n' = 0, and integrating with respect pz, we find 
the electron contribution into the dielectric function in 
the ultra-quantum limit 

n N 2e^pFi^c 

ee(k,w) = £oo -2 



{kanfl2 , {kaufl2 



(23) 



where p_f is the Fermi momentum in the magnetic field. 
In order to find this value, we assume that the electron 
concentration 

^!s\Xn{x ~ cpyjeH)'^ 

s 

is fixed by the defect concentration tiq. In the ultra- 
quantum limit, all the electrons are on the lowest level 
with n = 0. Integrating with respect py and Pz, we get 



and Eq. H23|l becomes 



(24) 



(25) 



ee(k, w) = - e 
(kaHf/2 ^ {kaHfl2 



where = Amru^e^ lm*eoo- Therefore Eq. H25(l coin- 
cides with Eq. lfTC|l . if l/r = and v\ = ^fiujc/2m* . We 
see, that the parameter A = {kvp /i^Jc)^ /5 ^ (k/H)^ in 
the semiclassical case, Eq. (|17|) . whereas the correspond- 
ing parameter in the ultra-quantum limit, Eq. 125(1 . is 
A' = {kaH)y2 - k^/H. 

The effect of collisions was not taken into account in 
Eq. (|25|l . It can be done with the help of the simple 
phenomenology similar to the paper 11]. But for the 
case of our interest l/r lu, when the coupled modes 
are observed, we can merely replace w by w -|- l/r in Eq. 
(ESI). 



V. THEORETICAL RAMAN SPECTRA 

The zeros of the dielectric function e(k, w), Eq. 
give the frequencies of the longitudinal modes. With the 
neglect of the electron and phonon F damping, the 
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FIG. 1: Frequencies of the phonon-plasmon coupled modes 
as functions of the carrier cyclotron frequency for the carrier 
plasma frequency ujp^ = 1.4 and cjlo = 1-1; all frequencies 
are in units ojto ■ The solid lines are obtain in solution of the 
cubic equation (see text), the dotted lines show the results 
ignoring the fc— dispersion. 



dielectric function, Eqs. Ijltil) and H25|l . is real and the 
frequencies of the corresponding vibrations are also real. 
In order to determine these frequencies, one has to solve 
a cubic equation. 

To a good approximation one can first omit the small 
term with k"^ and obtain two eigenfrequencies lo± 



FIG. 2; The same as in Fig. 1 but for Upe ~ 0.6. 



1 



UJ. 



TO 



' oj: 



pi 



■ UJ] 



peJ 



2 2 2i 2\2,^22 il/2 

±2 [(^c - ^TO - ^pi +^pe) + ^^pi^pe] ' 



(26) 



by solving the corresponding quadratic equation. The 
third frequency is found in the vicinity of the pole uj — 
2ujc of the omitted term. This approximation is shown 
in dotted lines in Figs. 1 and 2. 

The numerical results including the fc— dispersion are 
shown in solid lines. Now we must distinguish the semi- 
classical and ultra-quantum regimes. In the last case, 
the condition hujc > £f is fulfilled and the Fermi energy 
Ef depends on the magnetic field according to Eq. (Pl|l . 
Then the condition of ultra-quantum regime is rewritten 
as {2eH/hc)'^/^ > (27r)^rto for the fixed carrier concentra- 
tion uq. For instance, the ultra-quantum regime is real- 
ized for the concentration no = 10^^ cm~^ at H > 7.8 T 
and for no — 10^^ cm~'^ at H > 36 T. We suppose that 
the semiclassical condition, Eq. (|16() is obeyed in the in- 
terval of magnetic fields displaced in Fig. 1, whereas the 
ultra-quantum regime, Eq. H25I) is assumed in Fig. 2. 

With the help of Eqs. 10)) and l|HJl, we can plot the the- 
oretical Raman spectra for the large. Fig. 3 and low. Fig. 
4 electron concentrations. The Faust-Henry coefficient is 
taken C — —0.5 as is usually accepted in the absence of 
the magnetic field. We see always three peaks: the weak 
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FIG. 3: Calculated Raman spectra for the carrier plasma 
frequency and cyclotron frequency indicated in figure. All 
the frequencies are in units ljto; the phonon width F = 3 x 
10~^a;TO, the carrier relaxation rate = 2 x 10~^wto, the 
plasma frequency of ions ujpi = 0.4a;TO, and the momentum 
transfer kvp ~ O.Sujto- The cyclotron frequency is supposed 
to be in the semiclassical regime. 



resonance approximately at uj — 2ujc and two resonances 
corresponding uj± . Far from the crossing, the modes have 
mainly a character of the phonon or plasmon vibrations. 

From the splitting in Figs. 1 and 2, we can conclude 
that the resonance at uj = 2ujc interacts more intensively 
with the mode of the plasmon character (0;+ for the large 
concentration of carriers and uj^ for the low concentra- 
tion). 

The figure 3 (see, also. Fig. 1) clearly demonstrates the 
result of screening in the case of a relatively large con- 
centration of carriers: the frequency which equals 
the frequency of the longitudinal mode in the absence of 
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Raman shift (rel. un.) 

FIG. 4: The same as in Fig. 3 but for the carrier plasma 
frequency indicated in figure, = 10~^u)to, and for the 
ultra-quantum regime. 



at Lu = 2ujc., a phonon-plasmon crossover is seen in the 
curves for ujc = 0.75wto and I.Owto- As the magnetic 
field increases, the weaker plasmon peak, after passing 
through the phonon one, becomes more intense. 



VI. CONCLUSIONS 

In this work we investigated the spatial dispersion ef- 
fect in the Raman scattering on the phonon-plasmon 
coupled modes. In spite of the small value of corre- 
sponding parameters kvp/t^ (for the semiclassical mag- 
netic fields) and kan (for the ultra-quantum regime), 
this effect is observable in the high magnetic fields when 
the cyclotron resonances intersect the frequencies of the 
phonon-plasmon coupled modes. The magnetic field pro- 
vides an additional parameter which gives a possibility to 
compare the effect of the electron interaction with both 
the longitudinal and transverse optical modes in doped 
semiconductors. 



carriers (at the particular choice of parameters, Wlo = 
l.la;To)j turns out to be close to the frequency of the 
transverse mode Wto • The case of the relatively low car- 
rier concentration is illustrated by Figs. 2 and 4. It is 
seen in Fig. 4 that, in a weak field (wc = 0.25a;To), 
the frequency of the longitudinal mode is close (as it 
must be) to wlo = I-Iwto- Besides the weak resonance 
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